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We calculate total masses and radii of neutron stars (NS) for pure neutron matter and nuclear 
matter in /3-equilibrium. We apply a relativistic nuclear matter equation of state (EOS) derived 
from Dirac-Brueckner-Hartree-Fock (DBHF) calculations. We use realistic nucleon-nucleon (NN) 
interactions denned in the framework of the meson exchange potential models. Our results are 
compared with other theoretical predictions and recent observational data. Suggestions for further 
study are discussed. 



1. INTRODUCTION 

One of the most important and challenging problems 
in both theoretical and experimental nuclear physics is 
to understand the properties of matter under extreme 
conditions of density and pressure. The determination 
of the EOS (namely, the relationship between pressure 
and density) associated with such matter is a non-trivial 
problem which has attracted significant effort over the 
last few decades. Concerning terrestrial systems, detailed 
knowledge of the EOS is important, for instance, for un- 
derstanding heavy-ion collision dynamics. On the other 
hand, the EOS is crucial for determining the properties 
of one of the most exotic objects in the universe-neutron 
stars. Neutron star properties depend on the knowledge 
of the EOS over a wide range of densities - from the den- 
sity of iron at the stellar surface up to several times the 
density of normal nuclear matter in the core region of the 
star pj. 

Neutron stars are the smallest and densest stars known 
to exist. Like all stars neutron stars rotate-some at a rate 
of a few hundred revolutions per second. Such a fast ro- 
tating object will experience enormous centrifugal force 
which must be balanced by another force, gravity in this 
case, to prevent the star from falling apart. The balance 
of the two forces sets the lower limit on the stellar den- 
sity 0- Some neutron stars are in binary orbit with a 
companion and, in some cases, the application of orbital 
dynamics allows for an assessment of the masses. So far 
six neutron star binaries are known and all of them have 
masses in the range 1.36 ± 0.08M© Clearly, obser- 
vations of NS masses and radii impose important con- 
straints on the EOS of dense matter, as the latter consti- 
tutes the basic input quantity that enters the structure 
equations of a neutron star 0, 0, |a| . 

Although considerable progress has been made, the 
EOS of dense matter remains uncertain at densities above 
p n ~ 3 x 10 14 g/cm 3 |(|. For the reason stated above, 
this places uncertainty on the calculated NS properties. 
Determining the high density stellar matter EOS is a 
tremendous task. While at densities p ~ po, (po ~ 
0.17/m -3 is the normal nuclear density), matter con- 
sists mainly of nucleons and leptons, at higher densities 
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several other species of particles are expected to appear 
due to the rapid rise of the baryon chemical potentials 
with density [?|. Among these are strange baryons, such 
as the A, E, and H hyperons. In addition, pion-nucleon 
resonances may appear in stellar matter along with pion 
and kaon condensations. Moreover, at super high densi- 
ties nuclear matter is expected to undergo a transition to 
quark-gluon plasma. The value of the transition density 
may be obtained from QCD lattice calculations at finite 
baryon densities, but is presently still uncertain [3. 

After the first theoretical calculations of neutron star 
properties performed by Oppenheimer and Volkoff Q 
and independently by Tolman , many theoretical pre- 
dictions appeared in literature. Non-relativistic and 
relativistic approaches have been used. A number of 
early theoretical investigations on NS properties were 
done within the non-relativistic Skyrme framework. The 
reader is referred to the work of F. Douchin and P. 
Haensel in Ref. The conventional Brucckncr the- 

ory with a continuous choice for the single particle po- 
tential and three-body forces was applied in Ref. . 
Numerous predictions allow for hyperons together with 
nucleons and leptons in generalized /3-equilibrium to be 
included in the stellar interior. The Brueckner-Hartree- 
Fock (BHF) scheme was extended to include these con- 
tributions [T2I Il3| and applied in NS calculations 0. 
Since the Walecka model [lj] was proposed and applied 
to nuclear matter properties, the relativistic mean field 
approach has been widely used in determining NS total 
masses and radii Il5l The DBHF approach was used 
to compute NS properties in Refs. [(J [13 ■ 

In general, the EOS of stellar matter is considerably 
model dependent and, as a result, predictions of NS prop- 
erties are quite different from model to model. Among 
the sources of model dependence are [lj: (1) the many- 
body framework used for the determination of the EOS, 
(2) the model used for the bare NN interaction, (3) the 
kinds of hadrons/leptons included in the description of 
electrically neutral NS matter, (4) the treatment of pion 
and kaon condensations, (5) including the effects of fast 
rotation, (6) allowing for a phase transition from con- 
fined hadronic matter to deconfined quark matter. The 
possibility of such transition has attracted great interest 
over the last 30 years El E E3, H3- It is generally 
agreed upon that the high pressure found in the core 
of a neutron star creates an ideal physical environment 
for hadrons to transform into quark matter - a state of 
matter of practically infinite lifetime. Nevertheless, until 
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recently no empirical evidence had been reported. To a 
certain extent, this may signify that the occurrence (or 
not) of a transition to quark matter in the NS interior 
has only a minor impact on the overall static properties 
(range of possible masses, radii, limiting rotational pe- 
riods) of the star [lj. The situation, however, could be 
completely different for the timing structure of a pulsar 
(its braking behavior as it evolves) which is expected to 
diverge significantly from the normal pulsar behavior if 
the star converts a fraction of its core matter to pure 
quark matter 0. 

Recently we have been concerned with probing the be- 
havior of the isospin-asymmetric EOS [22|. In our work 
on neutron skins [23j , we have studied applications of the 
EOS at densities typical for normal nuclei. In Ref. 24] we 
discussed the one-body potentials for protons and neu- 
trons from DBHF calculations of neutron rich matter (in 
particular their dependence upon the degree of asym- 
metry in proton and neutron concentrations). More re- 
cently, we applied our microscopic approach in calcula- 
tions of effective in- medium NN cross sections j2j| . It is 
also important and timely to look into systems which are 
likely to constrain the behavior of the EOS at higher den- 
sities, where the largest model dependence is observed. 
This is the main purpose of the present paper, where we 
report predictions for masses and radii of static (non- 
rotating) neutron stars. Aside from minor adjustments, 
the framework is the one described in Ref. . We apply 
a relativistic (DBHF-basccl) EOS and perform our calcu- 
lations for two distinct cases: (1) pure neutron matter, 
and (2) baryon/lepton matter in /3-equilibrium. It is also 
our objective to provide an overview of the present status 
of both theory and experimental constraints. Together 
with systematic calculations based on our microscopic 
model, this broad outlook will help us gauge the quality 
of our tools and determine the importance of potentially 
missing mechanisms. 

This work is organized in the following way: after 
the introductory notes in this section, we briefly re- 
view our theoretical framework (section 2); the relativis- 
tic (DBHF-based) EOS for both neutron matter and (3- 
stable matter is discussed in section 3; our numerical re- 
sults are presented and discussed in section 4; we con- 
clude in section 5 with a short summary and suggestions 
for further studies. 



2. FORMALISM 

In order to provide the reader with a self-contained 
manuscript, in this section we outline briefly the formal- 
ism used to derive the isospin-asymmetric EOS. A more 
extensive discussion can be found in Ref. [2^ and refer- 
ences therein. 



2.1. Realistic nucleon-nucleon interactions 

The starting point of any microscopic calculation of 
nuclear structure or reactions is a realistic free-space NN 
interaction. A realistic and quantitative model for the 
nuclear force that has a reasonable basis in theory is the 



one-boson-exchange (OBE) model [2(|. In this frame- 
work bosons with masses below the nucleon mass are 
typically included. The model we apply in the present 
study consists of six bosons four of which are of major 
importance: (1) The pseudoscalar pion with a mass of 
about 138 MeV. It is the lightest meson and provides 
the long-range part of the NN potential and most of the 
tensor force; (2) The p vector meson - a 2ir P-wave res- 
onance with a mass of about 770 MeV and spin 1. Its 
major role is to cut down the pion tensor force at short 
range; (3) The u) vector meson - a 37r resonance of about 
783 MeV and spin 1. It creates a strong repulsive central 
force of short range and the short-ranged spin-orbit force; 
(4) The isoscalar-scalar a boson with a mass of about 550 
MeV. It provides the intermediate range attraction nec- 
essary for nuclear binding and can be understood as a 
simulation of the correlated S-wave 27r-exchange. 

The mesons are coupled to the nucleon through the fol- 
lowing meson-nucleon Lagrangians for the pseudoscalar 
(ps), scalar (s), and vector (v) fields, respectively: 

£ P v = -^-^yVW " 5 (l) 

TYlp S 

C s = g s H^ (s) (2) 



C v = -gvH 1 ^ ~ ffi^W^ ~ d^ ] ) (3) 

with ip the nucleon and (jffi the meson fields respectively 
(notation and conventions as in Ref. J27|). For isovector 
(isospin 1) mesons (such as p and 7r), is to be re- 
placed by t ■ <// Q ) , with t the usual Pauli matrices. 

From the above Lagrangians the OBE amplitudes can 
be derived (see, for instance, Ref. |2^|). The OBE poten- 
tial is defined as the sum of the OBE amplitudes of all 
exchanged mesons. 

In this work, we use the Bonn A, B, and C parameter- 
izations of the OBE potential formulated in the frame- 
work of the Thompson equation. The main difference 
between these three potentials is in the strength of the 
tensor force as reflected in the predicted D-state proba- 
bility of the deuteron, Po j2^. Bonn A has the weakest 
tensor force with P D = 4.5% (see Table 2 in Ref. [2SjV 
The Bonn B and C potentials predict 5.1% and 5.5%, 
respectively. It is well known that the strength of the 
tensor force is a crucial factor in determining the loca- 
tion of the nuclear matter saturation point on the Coester 
band H. 

2.2. Conventional Brueckner Theory 

We review here the main steps leading to the self- 
consistent calculation of the energy per particle in infi- 
nite nuclear matter, with or without isospin asymmetry, 
within the BHF approach. This is mainly to provide a 
baseline for the next subsection where we describe its 
relativistic extension (which is the method we actually 
apply). For a complete description of the BHF method 
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the reader is referred to Refs. |30|, |31|, |32| . For a review, 
see, for example, Ref. p^| . 

Nuclear matter is characterized by its total density or 
Fermi momentum Uf (and the degree of isospin asym- 
metry, in case of unequal proton and neutron densities). 
With fci and being the momenta of two nucleons 
with respect to nuclear matter, it is customary to formu- 
late the problem in terms of their relative momentum, 
K = |(ki — k^), and one- half of the center-of-mass mo- 
mentum, P = i(ki + k 2 ). (Clearly, k lj2 = P ± K.) 

The effective nucleon-nucleon interaction in nuclear 
matter is described in terms of the reaction matrix, or 
the Brueckner G-matrix, which satisfies the in-medium 
scattering equation (the Bethe-Goldstone equation) 

r d 3 K' 

G 4 ,(P;K,K ) - fy(K,Ko)-y — ^-(K,K') 
Q i3 -(K',P; fcjQG^- (P;K',K 



e*,(P,KO-(e*,)o(P,K ) 



(4) 



where ij—nn, pp, or np, and the asterisk signifies that 
medium effects are applied to these quantities. K , K, 
and K' are the initial, final, and the intermediate mo- 
menta, respectively. Vij is the two body OBE potential 
briefly discussed in the previous subsection. e*- is the en- 
ergy of the two-nucleon system, and (c*j)o is the starting 
energy. Thus, 



e* 3 (P,K) = e *(P,K) + e *(P,K), 



(5) 



with e* j the total energy of a single nucleon in nuclear 
matter. 

The G-matrix equation (Eq. (4)) is density dependent 
due to the presence of the Pauli projection operator, Q, 
which prevents scattering into occupied states and is de- 
fined as 



Q y (K,P,fc F ) 



1 if ki > k F and ki > k F 
otherwise. 



(6) 



Detailed expressions for the (angle-averaged) Pauli oper- 
ator can be found, for instance, in Ref. |23 for the case 
of unequal proton and neutron Fermi momenta. 

Equation (4) is density dependent also through the 
single-particle energy, 



Ti(p) + Ui(p), 



(7) 



with Tj (p) the kinetic energy and [/, (p) the potential en- 
ergy due to the interaction of the nucleon with all the 
others in the medium. We define 



Ui(p) =< p\Ui\p > = Re[ < pq\ G in \pq - qp > 

q<k™ 

E 



+ 



< Pq\G ip \pq - qp >] (8) 



with \p > and \q > single particle momentum, spin, and 
isospin states. Schematically, Eq. (8) has the form 



U n = U n 



Up — Upp 



u, 



n p 



(9) 



(10) 



for neutrons and protons, respectively. Together with the 
G-matrix equation, Eqs. (9-10) constitute a coupled, self- 
consistency problem, which is solved using the "effective 
mass approximation" |3l|. 

To demonstrate the self-consistency procedure, we take 
the non-relativistic single-particle energy and set 



Zrrii 1ml 



O.i . 



(11) 



which amounts to parametrizing the potential Uiijp) in 
terms of the nucleon effective mass m* and the constant 
Uq^ (i = n,p). Equation (11) implies 



U i(P) = o * 

2 m,mj 



-V + U . 



(12) 



Hence in the non-relativistic case the single-particle po- 
tential, Ui(p), has been fitted with a quadratic function 
of p. (The subscript "i" signifies that the parameters are 
different for neutrons and protons.) Starting with some 
initial values of m* and Uq^, the G-matrix equation is 
solved and a first solution for Ui{p) is obtained. These 
solutions are then parameterized in terms of a new set 
of constants, and the procedure is repeated until conver- 
gence is reached (that is, until differences between the pa- 
rameters from successive iterations are within the desired 
accuracy). Once the single-particle potential is available, 
we calculate the (average) energy per particle (to lowest 
order in the G-matrix). 

As is well known, nuclear matter calculations based on 
the BHF approach fail to predict correctly the satura- 
tion properties of nuclear matter. The typical trend is 
that the saturation density is too high for reasonable en- 
ergies. Different calculations based on the conventional 
BHF approach differ somewhat in their predictions of 
nuclear matter saturation properties, but all of them fail 
to predict simultaneously the correct saturation energy 
and density [26|. For this reason, it has become popular 
to implement non-relativistic calculations with contribu- 
tions from phcnomenological three-body forces. 

On the other hand, already more than two decades ago 
it was realized that the explicit treatment of the lower 
component of the Dirac spinor in the medium could pro- 
vide the missing saturation mechanism. This observation 
started what became known as the DBHF approach. We 
review some of the main points in the next section. 



2.3. The Dirac-Brueckner-Hartree-Fock approach 

For a detailed description see, for example, Refs. 0, 
l34l l35l | . The essential point of the DBHF approach is to 
use the Dirac equation for the single particle motion in 
nuclear matter 



mi - Ui{p))ui(p, s) = 0, 



(13) 



where the most gen eral Lorentz structure for Uiijp) is 
approximated as 



C<(p)«Efe,i(p)+7oE#,i(p). 



(14) 



Here Us,i is an attractive scalar field and Uy 4 is the time- 
like component of a repulsive vector field. The fields Us,i 
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and Uy i are in the order of several hundred MeV and 
strongly density dependent |3(j. With the definitions 
m*(p) =mi + UsAp) and (PiT =P°- Uy/p), the Dirac 
equation in nuclear matter can be written as 

(f t -m*) Ul (p,s) = 0. (15) 
The positive energy solution of Eq. (15) is given by 



Ui(p, s) 



E*(p)+m* 
2m* 



1/2 



1 



E* ,+m* 



Xs, 



(16) 



with xs a Pauli spinor, and E*(p) = (m* 2 + p 2 ) 1 / 2 . Wc 
notice that this is formally identical to the free-space 
spinor, but with m, replaced by m*. 

The single particle potential is then calculated from 
the operator given in Eq. (14) and properly normalized 
Dirac spinors. Namely, 

Tin* 

Ui{ P ) — -^7 < pWsM + i°KM\p > • 

The quantities Us,% and U Vi are momentum dependent, 
but to a good approximation can be taken as constant 



33, 35]. Hence, for the single particle potential, one can 
write 



U S ,i + Uy, 



(18) 



which, as in Eq. (12), again amounts to parameteriz- 
ing the single particle potential in terms of two constant 
quantities (for each type of nucleon). In analogy with 
the usual Hartree-Fock definition of the single particle 
potential (see Eq. (8)), we also write 



E*(p)E*(q) 



< Pq\9in\pq -qp> 



E 



m*m* 



< p E*(p)E*(q) 



< pq\gi P \pq - qp >] (19) 



where <?ij satisfies the relativistic scattering equation 



.g y (P,K,Ko)-^(K,K )-y — %( K,K) — e* 3 (P,K') - (e* J -) (P, K ) " 
Using the definitions = ^gij-^r and V*j = -^rVij-^t, Eq. (20) can be rewritten as 

g I3 (p;k,Ko) = ^(^ k o)- /S^'K-k v;:;;-^;- Z'r;,»vr - '-") 



d 3 K' g tJ (K^,P;fc f )G tJ (P;K',K ) 
(2^) 3 ij[ ' ' 4(P,K')-(c?,-)o(P,Ko) 



r 



which is formally identical to the non-relativistic G- 
matrix equation (Eq. (4)). We are then in a situation 
that is technically equivalent to the non-relativistic case, 
except for the Dirac structure of the single-particle en- 
ergy. Namely, Eqs. (19) and (21) (where the potential V*j 
is written in terms of in-medium spinors, Eq. (16)), must 
be solved self-consistently for G^ and the single particle 
potential, Ui (p) , with the help of the chosen parametriza- 
tion, Eq. (18). We then proceed to calculate the energy 
per particle. 

In summary, the most significant difference between 
the BHF and DBHF schemes is that in the latter case 
the nucleon wave function is obtained self-consistcntly, 
while in the BHF approach the free-space solution is used. 
This difference turns out to be an important one - as a 
result of the reduced nucleon mass in the medium, the 
lower component of the nucleon spinor is larger than in 
free space. This is well known to produce a density- 
dependent, repulsive-many body effect, with the result 
that the predicted saturation den sity and energy are con- 
sistent with the empirical values |33] ■ Physically, the de- 
scription of the nucleon in nuclear matter with a Dirac 
spinor as in Eq. (16) can be regarded as effectively tak- 
ing into account some many-body force contributions 



[37| (hence, the reduced need for inclusion of three-body 
forces in the relativistic scheme). 

2.4. Relativistic stellar structure equations 

In the present section we discuss the structure equa- 
tions of static (non-rotating) neutron stars. These are ob- 
jects of highly compressed matter so that the geometry of 
spacetime is changed considerably from that of flat space. 
Thus models of such stars need to be constructed within 
the framework of General Relativity combined with theo- 
ries of superdense matter. The connection between these 
two branches of physics is provided by Einstein's field 
equations 

G" u = 87rGT^(e,P(e)), (22) 

which couple the Einstein curvature tensor, G MI/ , to the 
energy-momentum density tensor, J 1 '"', of stellar matter. 
The tensor T^ v contains the equation of state in the form 
P(e) (pressure, P, as a function of the energy density, e), 
and G is the gravitational constant. 

Einstein's field equations are completely general and 
simple in appearance. However they are exceedingly 
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complicated to solve because of their non-linear character 
and because spacetime and matter act upon each other. 
There are only a few cases in which solutions can be found 
in closed form. One of the most important closed-form 
solutions is the Schwarzschild metric outside a spherical 
star. Another is the Kerr metric outside a rotating black 
hole. 

Starting from Einstein's field equations one can de- 
rive the structure equations of a static, spherically sym- 
metric, relativistic star. For an explicit treatment and 
derivation see for example the work of Fridolin Weber in 
Ref. [l| . These equations are generally known as Tolman- 
Oppenheimer-Volkoff (TOV) equations. In a system of 
units where c = G = 1, the TOV equations can be writ- 
ten as 



dP 
dr 



e(r)m(r) 



1 



P(r) 



e(r) 



47rr 3 P(r) 
m(r) 



2m(r) 



dm(r) 
dr 



<iTre(r)r 2 dr. 



(23) 



(24) 



To proceed to the solution of these equations, it is nec- 
essary to provide the EOS of stellar matter in the form 
P(e). Starting from some central energy density e c = e(0) 
at the center of the star (r = 0), and with the initial 
condition m(0) = 0, the above equations can be inte- 
grated outward until the pressure vanishes, signifying 
that the stellar edge is reached. Some care should be 
taken at r = since, as seen above, the TOV equa- 
tions are singular there. The point r = R where the 
pressure vanishes defines the the radius of the star and 
M = m(R) = An ^ e{r')r' 2 dr' its gravitational mass. 

For a given EOS, there is a unique relationship be- 
tween the stellar mass and the central density e c . Thus, 
for a particular EOS, there is a unique sequence of stars 
parameterized by the central density (or equivalently the 
central pressure P(0)). 

In this work we apply a standard fifth order Runge- 
Kutta numerical scheme [38) to integrate the TOV equa- 
tions, which are supplemented by the EOS in numerical 
form. In the next section, we discuss the features of the 
EOS applied in this study. 



3. EQUATION OF STATE 

We apply here the relativistic equation of state. For a 
detailed discussion of the isospin-dependent EOS we refer 
the reader to Ref. [z| • In the present study we calculate 
neutron star properties considering either pure neutron 
matter or asymmetric matter in /3-equilibrium. 



3.1. Neutron matter EOS 

The degree of isospin asymmetry is represented by the 
asymmetry parameter, a — — — — with p„ jP the neu- 
tron/proton densities and p = p n + p p the total density. 
(Clearly, a = for symmetric nuclear matter and a — 1 
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FIG. 1: (Color online) EOS for neutron (a = 1.0) and sym- 
metric matter (a = 0.0). 



for pure neutron matter.) The isopsin-dependent EOS 
can be written as 221 



e(kp, a) 



(1 + a)e n + (1 - a)e p 



(25) 



with e„ jP the average energy per neutron/proton and 
e(kp, a) the energy per particle as a function of the total 
Fermi momentum and the asymmetry parameter. The 
total Fermi momentum, kp, is related to the total den- 
sity, p, in the usual way 



Ikp 
3^" ' 



(26) 



We start with showing the isospin-dependent EOS for 
neutron and symmetric matter, see Fig. 1. In both cases, 
predictions are shown for DBHF and BHF calculations 
and all three NN potentials referred to in subsection 2.1. 
For the case of pure neutron matter (upper frame of 
Fig. 1), we see essentially no differences among the pre- 
dictions from the three potentials. As already discussed, 
a main source of differences among the three versions of 
the Bonn potential is in the strength of the tensor force, 
which is mostly reflected in the (T = 0) 3 Si — 3 D\ cou- 
pled states. In neutron matter (T = 1), however, this 
partial wave does not contribute and thus the model de- 
pendence is dramatically reduced j^. Consistent with 
these observations, in symmetric matter significant differ- 
ences are observed among the three potentials, see lower 
panel of Fig. 1. Notice that our EOSs are considerably 
more attractive than those shown in Ref. [2(j and also 
calculated with Bonn A, B, and C. For details on the dif- 
ferences between our respective calculations we refer the 
reader to Ref. |23| . 

For solving the TOV equations one needs to compute 
energy density and pressure, both of which are simply 
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FIG. 2: (Color online) Effective nucleon masses in neutron 
(a — 1.0) and symmetric (a = 0.0) matter as a function of 
density. For a — 0.0, the solid (black), dotted (dark blue), and 
dashed (light blue) curves correspond to Bonn A, B, and C, 
respectively. No differences can be seen among the predictions 
from the three potentials in the case of neutron matter. 



related to the EOS. As previously outlined, in our self- 
consistent calculation of the EOS we fit the single-particle 
potential with its "effective mass approximation" . Al- 
though chosen on reasonable theoretical grounds, the 
ansatz we use has a finite domain of validity (for which 
reliable convergenge of the self-consistent procedure can 
be obtained.) On the other hand, for NS calculations 
one needs to supply the relationship between energy den- 
sity and pressure over a very wide range of densities. 
The EOS at the lowest/highest densities have been de- 
rived with different methods, usually coupling available 
parametrizations of different EOSs Here, in the 

effort of keeping internal consistency as far as possible, 
we first obtain nucleon effective masses by interpolating 
the available predictions to the free-space value (to cover 
the small low-density part, approximately p < 0.03 fm -3 ) 
and extrapolating to the higher densities (approximately 
p > 0.6 fm -3 ). With the masses thus generated, we then 
proceed to the usual (microscopic) calculation of the en- 
ergy per particle everywhere in the needed range. (No- 
tice that the calculation of the G-matrix, and ultimately 
of the energy /particle, can proceed as long as effective 
masses are provided, as the momentum-independent part 
of the single-nucleon potential in Eq. (18) (or Eq. (12)) 
drops out from the energy denominator in Eq. (4).) 

The effective masses are shown in Fig. 2 for both sym- 
metric and neutron matter. We note the weaker den- 
sity dependence at the higher densities, a behavior that 
is already reflected in the self-consistent calculation and 
appears reasonable on physical grounds. The nucleon ef- 
fective masses originate from the (attractive) scalar po- 
tential in the Dirac equation, which, in turn, is sensitive 
mostly to the scalar meson a. As density increases, short- 
range contributions become dominant over intermediate- 
range attraction. This results into the observed lesser 
sensitivity of the effective masses to increasing density. 
(Further discussion on the high-density behavior of EOS 
will be presented in subsection 3.3.) 
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FIG. 3: (Color online) EOS for pure neutron matter. The 
upper panel shows the energy density and the lower panel the 
pressure as a function of the baryon number density. 

The energy density, e, is defined as 

e(p, a) = p[e(p, a) + m^ ee ] (27) 

where e(a, p) is the energy per nucleon and m^ L ' e is the 
nucleon rest mass. The pressure in nuclear matter is 
defined in terms of the energy per particle and the baryon 
number density as 

P(p,a)=p 2 ^. (28) 
dp 

In the upper frame of Fig. 3 we show the energy den- 
sity, e, as a function of the baryon number density, p, cal- 
culated with the Bonn A, B, and C potentials. The pres- 
sure, P, as a function of p is shown in the lower panel of 
Fig. 3. We also note that our predictions are reasonably 
consistent with most recent constraints on neutron and 
symmetric matter pressure obtained through medium en- 
ergy heavy ion collisions (Ijj. For the reasons discussed 
earlier, hardly any model dependence is seen in the pure 
T = 1 system. Accordingly, very small differences can be 
expected among the NS properties computed from pure 
neutron matter using the three potentials under consid- 
eration. (We reiterate that this is due to the relationship 
among these potentials in particular, and does not im- 
ply model independence in general.) Next we discuss the 
EOS for /3-stable matter. 



3.2. /^-stable matter EOS 

The density dependence of the symmetry energy deter- 
mines the proton fraction in /3-equilibrium, and, in turn, 
the cooling rate and neutrino emitting processes. 
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First, we consider a system of nucleons and electrons 
only. From the /3-st ability condition, the chemical poten- 
tial of the electrons is written as 



Pn Pp pe 



(29) 



Deriving the EOS of /3-stable matter requires knowing 
the corresponding proton fraction, Y p — —. The equilib- 
rium particle concentrations, Y^ = ^ (i = e, n, p), can 
be calculated via Eq. (29) combined with the charge neu- 
trality condition p e = p p (i.e. Y e = Y p ) and p = p n + p p 
(i.e. 1 = Y n + Y p ). To compute the various chemical 
potentials, 



Pi = 



dY t ' 



(30) 



we write etot as 



etot(p, Y p ,Y n ,Y e ) = e(p,Y p ) + —m p c 2 + —m n c 2 

P P 



Ee 

A ■ 

(31) 

Here e(p, Y p ) is the energy per nucleon and -f- is the con- 
tribution to the total energy from the electrons. Similarly 
to what is done in Ref. J41J, we assume a simple model 
for electrons and treat them as a gas of extremely rela- 
tivistic non-interacting fermions. That is, for electrons 
we set 



e{p 2 + m 2 e c 2 ) 1 ' 2 



pc, 



(32) 



which appears reasonable since the electron mass m e is 
small compared to its chemical potential at typical nu- 
clear density. Integrating over all momentum states (and 
summing over spin states), it is easy to show that the 
electron energy density (or energy per volume) becomes 



E e _ {k F f _ (3^ 2 pF e ) 4 / 3 
V 



4vr 2 



4tt 2 



(33) 



where kp is the electron Fermi momentum. 

In Ref. (2^ we have verified that the energy per nucleon 
can be written as 



e(k F ,a) = e(k F , 0) + e syrn (k F )a 2 , 



(34) 



with e(kp,0) the energy per particle for symmetric mat- 
ter and e sym the symmetry energy defined by 



_ 1 d 2 e(k Fl a) 



Q = 



With the definition a = 1 — 2Y p , Eq. (34) becomes 



e(p, Y p ) = e(p, Y P = \) + e sym (p)(l - 2Y p f 



Thus, etot is written as 



e to t(p, Y p , Y n , Y e ) = e(p, Y p = -) + 



+Y p m p c 2 + Y n m n c 2 + he 



4tt 2 



P 



(35) 



(36) 



2Yp) 2 



(37) 



where the last therm can be obtained trivially from 
Eq. (33) divided by the density. Knowing the total energy 



allows one to evaluate the chemical potentials, which, in 
turn, are inserted in Eq. (29) (combined with the con- 
straint Y p = Y e ) to give a simple algebraic equation for 

Yp 

-Ae sym {p)(l-2Yp) + {mp~m n )c 2 +hc{3n 2 p) 1 / 3 Y^ 3 = 0. 

(38) 

The solution of this equation is inserted in Eq. (36) to 
provide the EOS of nuclear matter in /3-cquilibrium with 
electrons. 

Just above the density of normal nuclear matter, the 
electron chemical potential, p e , exceeds the muon mass, 
m M , and the reaction n p + p~ becomes energetically 
allowed 0- Under these circumstances the /3-stability 
condition and charge-neutrality require 



and 



Pe — Pfit 



Pp = Pe + Pfi 



(39) 



(40) 



(and thus, Y p — Y e — Y^ = 0). Equation (31) needs to be 
modified to include the muon contribution to the total 
energy, 

etot(p, Y p , Y n , Y e , Yp) = e(p, Y p ) + Y p m p c 2 + Y n m n c 2 

E - E " (41) 



We will treat muons as a gas of non-interacting non- 
relativistic fermions |43|. a reasonable approximation 
given the relatively large muon mass (compared to its 
Fermi momentum at typical muon concentrations in (3- 
equilibrium). It can then be shown that the muon energy 
density is Q] 



2 n..W 3 



-T7- ~ P^rrip + 
V lvir z m b 



(^pp) 



(42) 



(in units where c = % = 1). Thus, Eq. (41) becomes 

1 

e to t(p, Y p ,Y n ,Y e , Yp) = e(p, Y p = -) + e sym (p){l - 2Y P Y 



+Y p in p c 2 + Y n m n c 2 + he- j-^- p 



(3^ e ) 4/3 1/3 
4tt 2 



+Y u m u c 



(he) 2 (3tt%) 5 / 3 2/3 



^ 2m u c 2 hit 2 9 ' 



(43) 



where the muon contribution to the total energy per par- 
ticle is obtained from Eq. (42) divided by p. As done in 
case of nucleons and electrons only, evaluating the chem- 
ical potentials, combined with Eqs. (39) and (40) leads 
to the following algebraic equations for the equilibrium 
particle concentrations (Yi) 

-4e sym (p)(l-2Yp) + (m p -m n )c 2 + hc(3Tr 2 p) 1 / 3 Y e 1 / 3 = 

(44) 



Vc 2 + (^) 2 ^P^P 2/3 Y 2/3 ~ M^ 2 P) 1/3 Y C 1/3 = 

(45) 

(46) 



Y p -Y e -Y il = 0. 
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FIG. 4: (Color online) EOS for neutron matter and 
baryon/lepton matter in (3 equilibrium. Predictions are ob- 
tained with the Bonn B potential. The upper panel shows the 
energy density and the lower panel the pressure as a function 
of the baryon number density. 



Once Y p is obtained from the above equations, it is in- 
serted in Eq. (36) to provide again the EOS of /3-stable 
matter. The energy density of baryon/lepton matter in 
/3-equilibrium is then used as the new input for the TOV 
equations. 

In Fig. 4 we compare EOSs for neutron matter and 
/3-stable matter (for both e~ and e~ plus \i~ cases), 
obtained with the Bonn B potential. We observe only 
very small differences among the various predictions, al- 
though the /3-stable EOSs are slightly "softer" than the 
one for pure neutron matter (at intermediate densities), 
a trend that is independent of the particular NN poten- 
tial being used. This is to be expected, since a system 
with some T=0 component is generally more attractive 
(mostly through the 3 Si partial wave) than a pure T=l 
system. This point will be looked at more closely in the 
next subsection. 

The symmetry energy and the proton fraction in /3- 
stable matter are displayed in Fig. 5. We observe that 
the symmetry energy (upper frame) tends to saturate 
and eventually decreases with density, a behavior qual- 
itatively similar to the one seen in the proton fraction 
(middle frame). From the lower panel of Fig. 5, we see 
that the addition of muons increases the proton fraction 
noticeably. In particular, the maximum proton fraction 
increases from 0.11 to about 0.13 (with Bonn B). We 
recall that a proton fraction larger than approximately 
^ would allow a critical cooling mechanism for the neu- 
tron star via the direct URCA processes ^4|. On the 
other hand, due to the behavior of the predicted sym- 
metry energy (the crucial player in the determination of 
the proton fraction), our proton concentrations remain 
very small, and thus the EOS is only mildly impacted 
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FIG. 5: (Color online) Symmetry energy, e sym , and proton 
fraction in /3-equilibrium, Yp. The upper panel shows the 
symmetry energy and the middle/lower panel the proton frac- 
tion as a function of density. 



by the inclusion of leptons. In Fig. 6 we show the differ- 
ent particle concentrations as a function of density. The 
presence of muons becomes possible only at higher den- 
sities, due to the energetic constraints mentioned earlier. 
Also, their concentration remains substantially below the 
electron fraction. 

In closing this section, we stress that the high-density 
behavior of the symmetry energy is very poorly con- 
strained and theoretically controversial, with different 
models often yielding dramatically different predictions. 
A more detailed discussion on this and related issues is 
presented next. 



3.3. High-density equation of state and the meson 
model: further discussion 



In this subsection, we take a more in-depth look at 
how the observed features of the EOS and the symme- 
try energy can be understood in relation to the physical 
components of our model. 

First, we observe that from Eq. (34) one can write 



e(k Fl a = 1) - e(k F , a = 0) = e sym (k F ) 



(47) 



which clearly displays the significance of the symmetry 
energy as the energy shift between neutron and symmet- 
ric matter. Thus, the behavior of e sym reflects the re- 
lationship between e(k F ,a — 1) and e(k F ,a = 0) as 
a function of density. What we observe in the present 
calculation is that the energy per particle in symmetric 
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FIG. 6: (Color online) Concentrations of various particle 
species (Yi = pi/ p) in /3-stable matter as a function of the 
baryon number density in units of po. 



matter becomes more and more repulsive with increasing 
density, eventually approaching the one of neutron mat- 
ter. Within the meson model, this can be understood in 
terms of the competing roles of the intermediate-range at- 
traction (provided by the scalar meson a and the iterated 
one-pion exchange) and the short-range repulsion (gen- 
erated by the vector meson to), as we describe next. Neu- 
tron matter is generally a more repulsive system because 
it lacks the contribution from T=0 waves, some of which 
are very attractive at normal nuclear densities (hence, 
for instance, the crucial role of the 3 <Si state in nuclear 
binding). On the other hand, as density increases, short- 
range repulsion becomes the dominant contribution to 
the NN interaction. Furthermore, in the DBHF calcula- 
tions the (repulsive) spin-orbit interaction is strongly en- 
hanced due to the relativistic effective mass in the OBE 
potential. Therefore, keeping in mind that only T=l 
states contribute to the energy of neutron matter while 
both isospin states contribute to the energy of nuclear 
matter, if major T=0 partial waves become increasingly 
repulsive at short distances, it is possible for the energy 
of symmetric matter to grow at a faster rate and eventu- 
ally approach the neutron matter EOS. This is just what 
we observe in our model, as can be seen from Fig. 7, 
where we compare the energy per particle in symmetric 
and neutron matter up to very high densities. (The ar- 
guments above also explain why the EOSs for neutron 
matter and /3-stable matter are very close to each other, 
cf. Fig. 4.) 

To explore this further as well as check internal con- 
sistency, we have performed some diagnostic tests. We 
find that the contributions to the average potential en- 
ergy of symmetric matter from only 1 Sq are -18.75 MeV, 
-22.75 MeV, and 5.554 MeV at fc F =1.4, 2.0, and 2.5 
fm _1 , respectively. Similarly including only the 3 Si and 
3 Di states, we find (for the same Fermi momenta) -18.74 
MeV, -10.04 MeV, and 18.16 MeV. A correct way to state 
the result of the above tests is to say that, in the presence 
of repulsive forces only (a scenario simulated by the pres- 
ence of high-density), nuclear matter would be a more 
repulsive system than neutron matter (for the same Uf)- 

At some critical density, signaled by the symmetry en- 
ergy turning negative, it would then be possible for a pure 
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FIG. 7: (Color online) Energy per nucleon for neutron (red 
curve) and symmetric (black curve) matter. 



neutron system to become more stable than symmetric 
matter, a phenomenon referred to as isospin separation 
instability 45] . Clearly, the value of such critical den- 
sity depends upon the relative degrees of attraction and 
repulsion in the particular model under consideration. 
(Figure 5 suggests that in the present model this would 
happen at densities well above ten times nuclear matter 
density.) 

It is of course possible that contributions not included 
in our model would soften the EOS in such a way as to 
alter the balance between the curves shown in Fig. 7 and, 
in turn, the high-density behavior of the symmetry en- 
ergy. This is precisely among the aspects we wish to learn 
about with this and future studies (namely, domain of va- 
lidity of the meson-theoretic picture, density dependence 
of the repulsive core, etc.). Systematic calculations 
based not on phenomenology but on a consistent the- 
oretical framework (in spite of its inherent limitations), 
together with stringent constraints, can help achieve that 
purpose. Our main conclusion at this point is that empir- 
ical constraints specifically on the high-density behavior 
of the symmetry energy would provide some clear and di- 
rect information on the short-range nature of the nuclear 
force. 



4. RESULTS FOR NEUTRON STAR TOTAL 
MASSES AND RADII 

In Fig. 8 we show the total masses and radii versus 
the stellar central density, p c . As expected in the light 
of our previous discussion, the NS properties computed 
from the EOS of neutron matter are essentially the same 
for all of the three potentials applied here. Consistent 
with the observations contained in subsection 3.3, the (3- 
stable matter EOS, being only slightly softer than the 
pure neutron matter EOS, yields similar predictions, al- 
though radii tend to be smaller at higher central densi- 
ties. In summary the neutron matter EOS from DBHF 
calculations yields a maximum mass M max ~ 2.25 Mq 
at radius R ~ 11.01 km and central density p c ~ 0.98 
/to -3 . The relativistic EOS for /3-stable matter yields 
maximum masses in the range M max w (2.238 — 2.241) 
M Q at radii R w (10.74 - 10.87) km and central densi- 
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ties p c w (0.998 - 1.013) /to" 3 . Including muons in the 
description of /3-stable matter does not alter significantly 
the predicted NS properties from the electrons-only case. 
Table 1 provides an overview of our predictions. Figure 
9 displays the mass-radius relation for the NS models 
shown in Fig. 8. 

Before we proceed with our discussion, we recall (see 
section 3) that an extrapolation is applied to the nucleon 
effective mass at those densities where a self-consistent 
solution for the single-particle potential is not available. 
Of course, every numerical extrapolation/interpolation 
involves uncertainties of some degree. To investigate to 
what extent the NS properties are sensitive to variations 
of to*, we have performed numerous tests all of which 
verify that the predicted NS masses and radii are not af- 
fected significantly by moderate variations of the nucleon 
effective mass. For instance, uncertainty of 20% in the 
high density effective mass results in approximately 2.9% 
variation in the NS maximum mass, 1.54% in the radius, 
and about 0.31% in the central density (e.g. Bonn B, 
DBHF, neutron matter: M max = 2.2457 ± O.O659M , 
R = 11.001 ± 0.17 km, p c = 0.979 ± 0.003/to' 3 ). Very 
similar conclusions apply for the case of /3-stable matter. 
(Sensitivity tests are performed in the following way: the 
highest density mass in Fig. 2 is varied by up to 20% and 
the self-consistcntly calculated values are then joined to 
it by interpolation. In this way the applied mass varia- 
tion increases gradually as we move to higher densities, 
where the uncertainty should be largest.) 

Our findings are generally consistent with previous 
ones obtained from similar theoretical frameworks [(J 0, 
E3- In Ref. @, for instance, both BHF and DBHF ap- 
proaches are used together with the Bonn A potential. 
Using the relativistic EOS, the authors predict a maxi- 
mum mass M max m 2.37 Mq with a radius R rs 12.2 km 
at a central density of p c « 0.8 /m~ 3 . We must keep 
in mind that the radius of a neutron star is mostly sen- 
sitive to the difference in pressure between neutron and 
symmetric matter |46) (the same mechanism that pushes 
neutrons out in the skin of a large nucleus). We also 
recall that, although our EOS's are rather repulsive (a 
feature generally shared by DBHF predictions), for the 
reasons discussed in subsection 3.3 the symmetry energy 
does not keep on growing with increasing density. As a 
consequence, the radii we predict are smaller than, for 
instance, those reported in Ref. [(J. From Fig. 9, we 



TABLE 1: NS maximum masses, radii, and central densities 
(from DBHF calculations). 

potential model composition M rnax {MQ) R(km) p c (fm~ 3 ) 
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FIG. 8: (Color online) NS total masses and radii. In all frames 
the solid (black) curve corresponds to the DBHF calculation 
for pure neutron matter, the short-dashed (red) curve to the 
DBHF calculation for /3-stable matter of nucleons and elec- 
trons only, and the long-dashed (blue) curve to the DBHF cal- 
culation for nucleons, electrons and muons in /3-equilibrium. 
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FIG. 9: (Color online) Mass-radius relation. The curves are 
labeled as in Fig. 8. 



observe that our predictions are compatible with recent 
findings from analysis of heavy-ion collision observables 
[46j which constrain the radius of a 1.4 Mq neutron star 
to be between 11.5 and 13.6 km. 

The ultimate upper limit of the NS mass can be de- 
duced theoretically on the basis 01 that (1) general rela- 
tivity is the correct theory of gravity, (2) the EOS of stel- 
lar matter is known below some matching density, and 
(3) the EOS of NS matter satisfies both (i) the causal- 
ity condition Qf- < c 2 and (ii) the microscopic stability 



condition > (known as the Le Chatelie's principle). 
Under the above assumptions the upper mass limit was 
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found to be 3.2 Mq. It is generally accepted that stars 
with masses above 3.2 M© collapse to black holes. On the 
other hand, there is a practical theoretical lower limit for 
the NS gravitational mass. The minimum possible mass 
is evaluated to be about (1.1 — 1.2) M© and follows from 
the minimum mass of a proto-neutron star. This is es- 
timated by examining a lepton-rich configuration with 
a low-entropy inner core of approximately 0.6 M© and 
a high-entropy envelope |4S|| . This argument is in gen- 
eral agreement with the theoretical results of supernova 
calculations in which the inner collapsing core material 
comprises at least 1.0 M© |49j . 

On the experimental side, the masses of neutron stars 
are mainly deduced by observations of NS binary sys- 
tems. Precise measurements of the masses of the binary 
pulsar PSR 1936+16 yield 1.344 M© and 1.444 ± 0.008 
Mq [H(j ■ Recently a few other binary pulsars have been 
observed all having masses in the narrow range 1.35 ±0.04 
Mq [53|. In addition, there are several X-ray binary 
masses that have been measured to be well above the 
average of 1.4 M©. The non-relativistic pulsar PSR 
J1012+5307 is believed to have a mass of approximately 
2.35 ± 0.85 Mq [13. Also the mass of the Vela X-l pul- 
sar has been deduced to be approximately 1.9 ± 0.2 M© 
Another object of the same type is Cygnus X-2 
with a mass of 1.78 ± 0.23 M© j54[. More recently, the 
mass of the binary millisecond pulsar PSR J0751+1807 
has been measured (by relativistic orbital decay) to be 
2.1±0.2M© [5f|, which makes this neutron star the most 
massive ever detected. 

Although accurate masses of several neutron stars are 
available, precise measurements of the radii are not yet 
available. (As mentioned earlier, terrestrial nuclear labo- 
ratory data are presently being used to obtain constraints 
on NS radii.) It has been shown that the causality condi- 
tion can be used [5|J to set the lower limit of the radius 
to about 4.5 km. In general, estimates of NS radii from 
observations have given a wide range of results. Perhaps 
the most reliable estimates can be done on the basis of 
thermal emission observations from neutron stars sur- 
faces which yield values of the so-called radiation radius 



R a 



R 



1 2GM \ 2 
1 ~ -R^-) 



(48) 



a quantity related to the red-shift of the star's luminosity 
and temperature ^| . The authors of Ref. [57j give values 
of i?oo < 9.5 km for a black-body source and R < 10.0 km 
for an object with a magnetized H atmosphere. The dis- 
covery of the quasi-periodic oscillations from X-ray emit- 
ting neutron stars in binary systems provides a possible 
way of constraining the NS masses and radii. Einstein's 
general relativity predicts the existence of a maximum 
orbital frequency which yields a mass of approximately 
1.78 Mq at a radius of about 8.86 km. Including cor- 
rections due to stellar rotation induces small changes in 
these numbers !5N . Recently, a radio pulsar spinning at 
716 Hz, the fastest spinning neutron star to date, has 
been discovered [60|. The authors of Ref. |6(j conclude 
that if the pulsar has a mass less than 2 Mq, then its 
radius is constrained by its spin rate to be < 16 km. 
Finally, the latest observations from EXO 0748-676 are 
reported to be consistent with the constraints 2.10±0.28 



Mq and 13.8 ± 1.8 km for the mass and the radius, which 
favor relatively stiff equations of state [61| . 

In view of the above survey of presently available ex- 
perimental constraints and/or estimates, we may con- 
clude that our relativistic predictions are reasonable, al- 
though our values for the maximum masses are on the 
high side of the presently accepted range. On the other 
hand, it must be kept in mind that mechanisms not in- 
cluded in the present calculation but generally agreed to 
take place in the star interior would further soften the 
EOS, thus decreasing the computed NS mass and radius. 
Among these mechanisms are pion and kaon condensa- 
tions, increase in the hyperon population due to rising 
chemical potentials with density, and transition to quark 
matter. All these phenomena tend to lower the energy 
per particle, and decrease the NS mass and radius while 
rising the stellar central density. Considering the effect of 
rotation in the NS properties calculation would increase 
the mass (by about 15%) and so this repulsive contribu- 
tion may be in part "cancelled" by any of the attractive 
mechanisms mentioned above. In summary, it seems that 
inclusion of additional degrees of freedom would, overall, 
move our DBHF predictions in the right direction. Our 
conclusions are summarized in the next section. 



5. CONCLUSIONS 

We have presented systematic calculations of NS lim- 
iting masses and radii using relativistic EOSs. We have 
considered the case of pure neutron matter and nuclear 
matter in /3-cquilibrium. The NS properties obtained 
in each case are very similar to each other, a behavior 
closely related to the predicted density dependence of 
the symmetry energy and the proton fraction in chemi- 
cal equilibrium. The present analysis helped us establish 
a clear correlation between the high-density behavior of 
the symmetry energy and the nature of the repulsive core 
as it manifests itself in specific partial waves. We have 
discussed this issue in considerable details and stressed 
the importance of reliable experimental information to 
set stringent constraints for theoretical predictions of the 
symmetry energy at high density. 

After overviewing presently available empirical infor- 
mation, we conclude that our DBHF-based predictions 
of NS maximum masses are somewhat high, indicating 
that repulsion in the high-density EOS should be reduced 
(again, keeping in mind that the difference in pressure 
between neutron and symmetric matter will primarily 
control the radii). 

It is likely that other mechanisms, not included in the 
present hadronic model, may take place in the stellar 
interior which are likely to soften the EOS. Pions and 
kaons may be likely to condensate in the interior of neu- 
tron stars [5j| |U • It is generally agreed upon that 
pion/kaon populations increase the proton fraction and 
might cause a rapid cooling via the direct URCA process. 
The stellar matter EOS could soften considerably due to 
pion/kaon condensations, since these mechanisms tend 
to lower the symmetry energy. In addition, at higher 
densities different species of hyperons appear in the NS 
composition which also results in a softening of the stellar 
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matter EOS. The expected transition to quark matter at 
higher densities would lower the energy and soften EOS. 
Therefore, the present set of results give us confidence 
in our relativistic microscopic approach as a reasonable 
baseline from which to proceed. 

Another important issue is the effect magnetization 
may have on the EOS, namely, how the energy per parti- 
cle changes if matter in the stellar interior becomes spin- 
polarized. As a next step towards a better and more 
complete understanding of the physics of neutron stars, 
we are presently extending our framework to include the 
description of spin-polarized neutron/nuclear matter. 

Finally, the observations of compact stars will be 
greatly improved in the future by the Square Kilometer 
Array (SKA). The SKA is an internationally sponsored 



[1] F. Weber, Pulsars as Astrophysical Laboratories for Nu- 
clear and Particle Physics, IOP Publishing Ltd, Bristol 
and Philadelphia (1999). 

[2] N. K. Glendenning, Compact Stars, Nuclear Physics, 
Particle Physics, and General Relativity, Springer- Verlag 
New York, Inc. (1997). 

[3] S. E. Thorsett and D. Chakrabarty, Ap. J. 512 (1999) 
288. 

[4] P. Wang, S. Lawley, D. B. Leinweber, A. W. Thomas and 
A. G. Williams, Phys. Rev. C 72, 045801 (2005). 

[5] S. L. Shapiro and S. A. Feukolsky, Black Holes, White 
Dwarfs, and Neutron Stars, The Physics of the Compact 
Objects, John Wiley & Sons, Inc., New York (1983). 

[6] G. Bao, L. Engvik, M. Hjorth- Jensen, E. Osnes, E. 
0stgaard, Nucl. Phys. A 575 (1994) 707-732. 

[7] M. Baldo, G. F. Burgio and H.-J. Schulze, Phys. Rev. C 
61 55801 (2000). 

[8] J. Oppenheimer and G. Volkoff, Phys. Rev. 55, 374 
(1939). 

[9] R. C. Tolman, 1934, Proc. Nat. Acad. Sci. USA 20 
(1934); Phys. Rev. 55 (1939) 364. 
[10] F. Douchin and P. Haensel, Astron. and Astrophys. 380, 
151-167 (2001). 

[11] M. Baldo, G. F. Burgio, H. Q. Song, and F. Weber, Proc. 
Int. Workshop XXVI on Gross Properties of Nuclei and 
Nuclear Excitations, Hirshegg, Austria, January 11-17, 
1998. 

[12] M. Baldo, G. F. Burgio, and H. -J. Schulze, Phys. Rev. 

C 58, 3688 (1998). 
[13] H.-J. Schulze, A. Lejeune, J. Cugnon, M. Baldo, and U. 

Lombardo, Phys. Lett. B 355 21 (1995); Phys. Rev. C 

57, 704 (1998). 

[14] B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16 
(1986) 1. 

[15] S. F. Ban, J. Li, S. Q. Zhang, H. Y. Jia, J. P. Sang, and 

J. Meng, Phys. Rev. C 69 (2004). 
[16] S. Lawley, W. Bentz and A. W. Thomas, Nucl. Phys. 

Proc. Suppl. 141, 29 (2005). 
[17] K. Sumiyoshi, K. Oyamatsu, and H. Toki, Nucl. Phys. 

A595, 327-345 (1995). 
[18] T. Endo, T. Maruyam a, S. Chiba and T. Tatsumi, 

arXiv:hep-ph/0502216 
[19] G. F. Burgio, Nucl. Phys. A 749, 337 (2005). 
[20] C. F. Burgio, M. Baldo, P. K. Sahu, and H. -J. Schulze, 

Phys. Rev. C 66, 025802 (2002). 
[21] S. Pal, M. Hanauske, I. Zakout, H. Stocker, and W. 

Greiner, Phys. Rev. C 60, 015802 (1999). 
[22] D. Alonso and F. Sammarruca, Phys. Rev. C 67, 054301 



project with the goal to construct a radio-telescope with a 
total receiving surface of one million square meters. The 
SKA is a facility with a potential to detect from 10,000 to 
20,000 new pulsars, more than 1,000 millisecond pulsars 
and at least 100 compact relativistic binaries |64|. These 
future endeavors will provide a tool capable of probing 
the NS matter EOS at the extreme limits. 



Acknowledgements 

The authors acknowledge financial support from the U.S. 
Department of Energy under grant number DE-FG02- 
03ER41270. 



(2003). 

[23] D. Alonso and F. Sammarruca, Phys. Rev. C 68, 054305 

(2003) . 

[24] F. Sammarruca, W. Barredo, and P. Krastev, Phys. Rev. 

C 71, 064306 (2005). 
[25] F. Sammarruca and P. Krastev, Phys. Rev. C 73, 014001 

(2006). 

[26] R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989). 

[27] J. D. Bjorken and S. D. Drell, Relativistic Quantum Me- 
chanics, McGraw-Hill, New York (1964). 

[28] R. Brockmann and R. Machleidt, The Dirac Brueckner 
Approach, in: Int. Rev. Nucl. Phys., Vol.8, Nuclear Meth- 
ods and the Nuclear Equation of State, M. Baldo, ed. 
(World Scientific, Singapore, 1999) Chapter 2. 

[29] F. Coester, S. Cohen, B. D. Day, and C. M. Vincent, 
Phys. Rev. C 1, 769 (1970). 

[30] H. A. Bethe, Ann. Rev. Nucl. Sci. 21, 93 (1971). 

[31] M. I. Haftel and F. Tabakin, Nucl. Phys. A158, 1 (1970). 

[32] D. W. L. Sprung, Adv. Nucl. Phys. 5, 225 (1972). 

[33] R. Brockmann and R. Machleidt, Phys. Rev. C 42, 1965 
(1990). 

[34] C. J. Horowitz and B. D. Serot, Nucl. Phys. A464, 613 
(1987). 

[35] B. ter Haar and R. Malfliet, Phys. Rep. 149, 207 (1987). 

[36] R. Brockmann, RIKEN Rev. 26 (2000). 

[37] G.E. Brown, W. Weise, G. Baym, and J. Speth, Com- 
ments Nucl. Part. Phys. 17, 39 (1987). 

[38] W. H. Press, B. P. Flannery, S. A. Teukolsky, W. T. Vet- 
terling, Numerical Recipes in Fortran: The Art of Scien- 
tific Computing, Cambridge University Press (1992). 

[39] G. Q. Li, R. Machleidt, and R. Brockmann, Phys. Rev. 
C 46, 2782 (19 92). 

[40] P. Danielewicz, nucl-th/0512009 

[41] I. Bombaci and U. Lombardo, Phys. Rev. C 44, 1892- 
1900 (1991). 

[42] R. B. Wiringa, V. Fiks, and A. Fabrocini, Phys. Rev. C 

38, 1010-1037 (1988). 
[43] C. Maieron, M. Baldo, G.F. Burgio, and H.-J. Schulze, 

Phys. Rev. D 70, 043010 (2004). 
[44] J. M. Lattimer and M. Prakash, Science 304, 536-542 

(2004) . 

[45] B.-A. Li, Nucl. Phys. A 708 (2002) 365-390. 
[46] B.-A. Li and A. Steiner, jrmcl-th/05110641 
[47] C. E. Rhoades and R. Ruftini, Phys. Rev. Lett. 32, 324 
(1974). 

[48] J. -O. Goussard, P. Haensel, and J. L. Zdunik, Astron. 

and Astrophys. 330, 1005 (1998). 
[49] J. M. Lattimer and M. Prakash, Astrophys. J. 500, 426- 



13 



442 (2001). 

[50] M. Weisberg and J. H. Taylor, Phys. Rev. Lett. 52, 1348 
(1984). 

[51] S. E. Thorsett and D. Chakrabarty, Astrophys. J. 512, 
228 (1999). 

[52] M. H. van Kerkwijk, P. Bergeron and S. R. Kulkarni, 

Astrophys. J. 467, L89-L92 (1996). 
[53] H. Quaintrell el al, Astron. Astrophys. 401, 313 (2003). 
[54] A. R. King and H. Ritter, Mon. Not. R. Astron. Soc. 

309, 253-260 (1999). 
[55] D. J. Nice, E. M. Splaver, I. H. Stairs, O. Lohmer, A. 

Jessner, M. Kramer, and J. M. Cordes, Astrophys. J. 

634, 1242-1249 (2005). 
[56] J. M. Lattimer, M. Prakash, D. Masak, and A. Yahil, 



Astrophys. J. 355, 241 (1990). 
[57] A. Golden and A. Shearer, Astron. and Asprophys. 342, 
L5 (1999). 

[58] D. Psaltis el al, Astrophys. J. 501, L95 (1998). 

[59] K. Takahashi, Prog. Theor. Phys. 108, 689 (2002). 

[60] J. W. T. Hessels, S. M. Ransom, I. H. Stairs, 

P. C. C. Freire, V. M. Kaspi and F. Camilo, Science 311, 

1901-1904 (2006). 
[61] F. Ozel, Nature 441, 1115-1117 (2006). 
[62] T. Muto, Nucl. Phys. A 754, 350 (2005). 
[63] S. Kubis and M. Kutschera, Nucl. Phys. A 720, 189 

(2003). 

[64] M. Kramer, |ar"Xiv:ast ro-ph/0306456 



